Production of the LiBeB isotopes with 
time-dependent GCR fluxes 
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Abstract. 

D ■ It is well known that no time-independent galactic cosmic ray (GCR) flux model can 

| account for the observed linear relationship of Be and B abundances with metallicity. In this 

■ contribution we solve the time-dependent CR propagation equation using simple analytic 

CN| \ time dependences for the CR spectrum at the source, the galactic CR escape-length and gas 

density evolution. The CR fluxes so obtained are used to study the early evolution of the 
LiBeB abundances in the Galaxy. 

£q ; 1 A time-dependent model for the GCR propagation 



Recent observations of unevolved halo stars have confirmed the surprising linear re- 
lationship of the Be and B abundances with [Fe/H] (Q, [||). This primary behaviour 
5o ! °f Be and B cannot be understood on the basis of the standard model (SM) of LiBeB 

Q\ ■ production in the interstellar medium (ISM) from GCR spallation. This is because in 

the SM the LiBeB production rates are proportional to the global metallicity of the 
O ISM, so that a secondary behaviour with [Fe/H] is expected. However, the SM as- 

' \ sumes a GCR flux constant in time and equal to that currently oberved at the Earth. 

^ ■ There are many reasons to believe that the GCR flux was different in early epochs. 

\ In this contribution we solve the time-dependent propagation equation of CRs in the 

0$ • galaxy in a simple and parametrized way to test whether the SM can be compatible 

with observations when the hypothesis of the CR flux constancy in time is relaxed. 
We consider a general time-dependent equation for the GCR propagation: 



dF(E,t) F(E,t) d(w(E,t)F(E,t)) 



; dt A(E, t) dE 



q(E,t) (1) 



where F(E,t) is the number density of GCR particles per energy interval as a func- 
tion of energy per nucleon and time, A(E, t) is the loss time-scale including escape, 
A e {E,t), and inelastic nuclear interaction, A n (E,t). w(E,t) represents the ionization 
energy losses in the ISM and finally, q{E, t) is the source spectrum of GCRs. We 
make the following assumptions: 

1) The energy and time dependences of A e .„, w and q are separable: 

A e , n (E,t) = Al> n (E)A e 2 > n {t);w(E,t) = w 1 (E)w 2 (t);q{E,t) = qi(E)q 2 (t) (2) 

2) We assume that GCRs are produced by supernovae (SN) so that q 2 (t) — SNR{t), 
the SN rate, taken from a chemical evolution model which reproduces the current 
SN rates in the galaxy and the main characteristics of the solar neighbourhood. The 
source energy spectrum is: 

rv\ (E + E ) 

qi{E) = qo [E(E + 2E )]^ (3) 



with Eq = 938 MeV. The constant go is chosen in such a way that we can reproduce 
the observed current GCR flux at the Earth. 
3) A(E, t) can be split into two parts 
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= 1 (4) 

A(E,t) A e (E,t) A n (E,t) 1 ' 

where A e represents the escape time-scale from the galaxy, and A n is the loss time- 
scale against inelastic nuclear interaction. The energy dependence of A e is taken from 
[Q and that of A„ from || assuming a gas density of 1 cm~ 3 . 

4) The energy losses due to ionization in the ISM have been considered according 
to |[. 

5) The time-dependence of the ionization losses and of the inelastic nuclear inter- 
actions can be plausibly caused by temporal changes in the gas density. The time- 
dependence of the escape time-scale is more difficult to understand. Here we treat 
it as a free parameter. For both we assume, for the sake of simplicity, exponential 
functions: 

A|(i) = exp(-t/r esc ); A%{t) = exp(-t/T gas );w 2 (t) = exp(-t/r gas ) (5) 

where r gas is the characteristic e-folding time of gas exhaustion in the galaxy and 
r esc is the characteristic time-scale for variations of A e . We assume that A(E,t) was 
higher in the past, during the halo collapse (t< 2 Gyr). 

2 The Solution 

Eq. (1) can be solved using the method of characteristics. The solution can be 
simplified considerably if T gas ,T esc 3> 10 7 yr. Plausible values for T gas and r esc are 
longer than 10 7 yr, so in that case the solution is 




F ( E ^) = ... 7~Er\77. n\ I dE Qqi {E )q 2 (t~T loss (E ,t)) 

(6) 
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dE' 



Wi(E')Al(E')A%(t)Al(E') I 
where Ti oss (E' ,t) is the ionization energy loss time-scale calculated as: 



T loss (E',t)= [ dE" ! (7) 
Je wi(E")w 2 {t) 

Then, the rate of accumulation of LiBeB in the ISM is calculated in the usual way 
(see [||). We have considered the following cases: 

Case 1: The energy dependence of the GCR source spectrum is given by Eq. (3) 
with q 2 (t) equal to the SN rate evolution. The remaining parameters are constant in 
time. 

Case 2: Case 1 plus A 2 (t) from Eq. (5) with T esc — 0.5 Gyr during the first 2 
Gyr of galactic evolution and r esc = oo afterwards (i.e no time variation of A). This 
means that the CR escape time-scale at t = was a factor ~ 500 higher than now. 



log [ E(MeV/n) ] [Fe/H] 



-8.5 




[Fe/H] [Fe/H] 



Figure 1: Calculated GCR flux (upper left panel) and evolution of the LiBeB abundances 
in Case 1 (dashed line) compared with the predictions of the standard model (continuous 
line) . 



3 Results and conclusions 

Figures 1 and 2 show the calculated time-dependent GCR fluxes as well as the corre- 
sponding LiBeB evolution with [Fe/H]. The results are: 

Case 1: The flux is higher in the past and decreases slowly with time in the whole 
energy range. This higher flux implies a larger LiBeB production in the past than 
the SM model, although, this is still not enough to reproduce the observations. 

Case 2: The flux in the past is also higher and the spectrum is flatter at energies 
E> 100 MeV/n. In consequence the LiBeB production is even higher at early epochs 
than in Case 1. Note that the predicted Be evolution is nearly in agreement with the 
observations. 

On the other hand, the effect of considering the temporal evolution of the gas 
content in the galaxy is not very important for the LiBeB production. In this case 
the flux and spectrum are modified in a complex way depending on the value of T gas 
and the interplay between A e and A„. Our results show that variations of the GCR 
flux in time have to be considered to study the early production of LiBeB by spallation 
in the ISM. A higher flux and a flatter CR spectrum in the past are able to improve 
the agreement between the predictions of the (modified) SM with observations. This 
is evident for Be and B. For Li, however, another important source of production is 
still required. A higher flux in the past can be achieved by taking into account the 
temporal evolution of the GCR sources, which is naturally explained by the temporal 
variation of the SN rate, or by assuming that the escape time-scale A esc was higher in 
the past. The later assumption can improve the agreement between data and model 
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Figure 2: Calculated GCR flux (upper left panel) and evolution of the LiBeB abundances 
in Case 2 (dashed line) compared with the predictions of the standard model (continuous 
line) . 



considerably. However, it still requires a more solid theoretical justification (0). 
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